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The spectral properties of a variety of improved staggered operators are studied in quenched QCD. 
The systematic dependence of the infrared eigenvalue spectrum on i) improvement in the staggered 
operator, ii) improvement in the gauge field action, iii) lattice spacing and iv) lattice volume, is 
analyzed. It is observed that eigenmodes with small eigenvalues and large chirality appear as the 
level of improvement increases or as one approaches the continuum limit. These eigenmodes can be 
identified as the "zero modes" which contribute to the chirality associated, via the index theorem, 
with the topology of the background gauge field. This gives evidence that staggered fermions are 
sensitive to gauge field topology. After successfully identifying these would-be chiral zero modes, 
the distribution of the remaining non-chiral modes is compared with the predictions of Random 
Matrix Theory in different topological sectors. Satisfactory agreement is obtained. 



I. INTRODUCTION 

It is well known that spectral flow of the massless Dirac 
operator [D] is related to the topology of the background 
gauge fields. When massless fermions are coupled to 
Yang-Mills fields with non-trivial topology, one or more 
eigenvalues of the corresponding Dirac operator necessar- 
ily vanish. This is formally given by the Atiyah-Singer 
index theorem Q 



Q = index(_D) = — n^j. 



(1) 



where Q is the topological charge of the gauge field and 
nLinR) are the numbers of left(right) zero eigenmodes of 
the Dirac operator. 

On the other hand, until recently it was believed that 
staggered fermions do not feel gauge field topology 2] 
because of the lack of zero eigenmode of the operator at 
finite lattice spacing. This unpleasant feature of stag- 
gered fermions was further revealed in comparisons of 
the infrared eigenvalue spectrum obtained in simulations 
with the predictions of Random Matrix Theory (RMT) 
0, 0, I3- The eigenvalue spectrum in all topological 
charge sectors was found to be consistent with thepre- 
diction of RMT for topological charge equal to zero |j, I3 ■ 

In the past year this problem has been revisited and the 
situation has changed. Recent calculations 0, 13, show 
that with staggered operator improvement for staggered 
fermions the correct response to QCD topologj^governed 
by the index theorem, can be obtained.^ In [3| we pre- 
sented some of our results in a preliminary form. Here 
we give our final results along with the study of some 
systematic effects not included in our earlier work. 



^ Earlier Farchioni et al. Q suggested that in the Schwinger model 
topology may be visible in the spectrum of the staggered Dirac 
operator. 



The failure of unimproved staggered fermion actions 
to show the proper topological features is due to lat- 
tice artifacts, in particular flavour-ch ang ing interactions. 
This was already suggested in Ref. The staggered 

quark action describes four quark flavours in the con- 
tinuum limit and the eigenvalue spectrum has a 4-fold 
degeneracy in this limit. At finite lattice spacing, the 
flavour-changing interactions break the fiavour symme- 
try and the degeneracy is lifted. Consequently, staggered 
fermions do not have exact zero modes at finite lattice 
spacing because the continuum chiral modes (if there are 
any) are scattered on the lattice. One thus expects stag- 
gered fermions to show better topological properties if 
flavour-changing effects can be suppressed. 

In addition to lattice artifacts, the spectral density in 
the infrared limit also depends sensitively on the volume 
of the lattice. This can be seen from the fact that, be- 
cause the topological susceptibility deflned by {Q'^)/V 
is independent of the volume V , the ensemble average 
(Q^) scales approximately with V . For example, pre- 
vious studies [ill . [T2i] on the spectral properties of the 
overlap operator have revealed this strong flnite volume 
effect (in contrast with the staggered Dirac operator, the 
overlap fermion operator 14] has exact chiral symmetry 
on the lattice and has true zero modes even at finite 
lattice spacing, see section llIl|) . In particular, it has been 
shown that satisfactory agreement with the predictions 
of RMT can be achieved only when the physical volume 
of the lattice exceeds (1.2/m)^. On smaller volumes the 
agreement is less appealing and the sensitivity to topol- 
ogy disappears. 

Both of these issues have been examined in this project 
by studying the spectral flow in the infrared limit for a 
variety of improved staggered operators. The improved 
staggered operators are designed to suppress flavour- 
changing interactions by smoothing out the quark-gluon 
interaction vertex |l6j | . Four different sources of system- 
atic effects have been considered: 
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• improvement in the staggered operator 

• improvement in the gauge field action 

• lattice spacing 

• lattice volume. 

Simulations are done both with the standard Wilson 
gluon action and the 0{a^as) accurate tadpole improved 
gluon action 0| for a range of lattice spacings and vol- 
umes. It is observed that eigenmodes with small eigenval- 
ues and large chirality appear as the level of improvement 
increases 0, |1| • These small eigenmodes can be identi- 
fied as the chiral zero modes associated with the topology 
of the background gauge fields. Our results also show the 
expected strong dependence on volume. Only a few "zero 
modes" are observed when the physical volume is small 
and the number increases as the volume increases. 

The paper is organized as follows. The improved stag- 
gered operators and their properties will be reviewed in 
the next section. A similar review on the overlap operator 
is given in section Hill Details of the simulations are dis- 
cussed in section Hvl and results will be presented in sec- 
tion After successfully identifying the "zero modes" , 
distribution of the remaining non-chiral modes is com- 
pared with the predictions of RMT in Section IVll 

II. IMPROVED STAGGERED OPERATORS 

The unimproved staggered Dirac operator is 

(2) 

where ?7^(a;) = - jg the standard staggered 

fermion phase. A variety of improvement schemes are 
considered here. The basic improved version is the C(a^) 
tree-level improved Asq operator IJ^ , which includes an 
additional 3-link Naik term and replaces the gauge field 
in ^ with Fat7 effective links (sum of the original link 
and the nearby paths, up to 7-link staples). Tadpole im- 
provement of the Asq operator gives the 0{a^as) accu- 
rate Asqtad operator. This is the staggered fermion op- 
erator which has been used in recent dynamical fermion 
simulations [T^ . 

It has been found that there are alternate improve- 
ment schemes that can suppress fiavour-symmetry break- 
ing lattice artifacts even more than Asq improvement. 
One such scheme is based on unitarized Fat7 (UFatT) 
links. The usefulness of UFat7 smearing was first dis- 
cussed by Lee |23| and a staggered action which iter- 
ates this fattening procedure giving the improved Asq 
operators (UFatT)" x Asq was studied by FoUana et al. 
I2H . The HYP-improved operators (HYP)" are con- 
structed from unitarized fat links also but only those links 
within the hypercubes containing the original link are in- 
cluded in the fattening process. Readers are referred to 
the respective papers for further details. 



The staggered Dirac operator is anti-Hermitian and 
since {D',e{x)} = where e{x) = (_l)^i+ - +^4^ the 
eigenvalues of D'^ are pure imaginary and come in pairs 
±iA. The commutation relation {-D*, e(a^)} = is a rem- 
nant of the continuum global chiral symmetry on the lat- 
tice. In addition, it is well known that {D^Y connects 
only even-even ((-D|e)^) oi" odd-odd {{Dl^Y) sites on the 
lattice. In this project, we choose to compute the eigen- 
values of {D^^Y in favour of D'' since it is Hermitian and 
has real eigenvalues A^. Note that the extra doubling of 
modes induced by squaring is canceled by working only 
on even sites of the lattice. To be precise, we compute 
the lowest 40 eigenvalues of {Dl^Y which correspond to 
the lowest 40 positive (imaginary) eigenvalues of D" . 

Before leaving this section it is useful to understand 
how the low-lying non-chiral modes scale with the var- 
ious systematics. This explains the scales used in the 
spectral plots given later in section In Ref . ^3 j 
Kalkreuter computed the complete spectrum of the unim- 
proved staggered Dirac operator on 12^ lattices at several 
lattice spacings. Results show an accumulation of small 
eigenvalues with increasing lattice spacing. In addition, 
the number of eigenvalues increases when the volume of 
the lattice increases. Consequently, the magnitude of 
the non-zero eigenvalues in the infrared limit decreases 
with increasing lattice spacing and volume. On the other 
hand, when lattice spacing decreases, i.e., one goes closer 
to the continuum limit, an increase in the magnitude of 
the low-lying non-chiral modes is expected. The same 
should happen as we come closer to the continuum by 
improving the action. These scaling properties will be 
refiected in the choice of axes for the spectral graphs. 



III. THE OVERLAP OPERATOR 

The massless overlap operator is given by |0| 

Z?° = l+75e(i?^"), (3) 
where e{H^) is the matrix sign function 



(4) 



and H^' — 75/?"' is related to the standard Wilson Dirac 
operator 



(5) 



The hopping parameter k is related to the bare mass m 
by 



1 



-2m + 8 



(6) 
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Action 


2 

a a 


a (fm) 


Volume V 


5.85 


Wilson 


0.0748(7) 


0.123 




6.0 


Wilson 


0.0478(5) 


0.0981 




6.2 


Wilson 


0.0259(5) 


0.0722 


16'' 


8.26 


improved 


0.0724(5) 


0.121 


8^ lo^ l2^ 16-' 


8.62 


improved 


0.0456(4) 


0.0958 




9.18 


improved 


0.0246(5) 


0.0704 


16" 



TABLE I: Simulation parameters and measured string ten- 
sions a^a. The couplings are chosen such that lattice spacings 
agree between simulations using the Wilson gluon action and 
the improved gauge field action. 



A detail description of how k should be chosen is given 
in ■ We use k = 0.21 or to = 1.62, which has been 
shown to be appropriate for our study of topology |3 • 

The overlap operator ^ satisfies the Ginsparg- Wilson 
relation IHl 



(7) 



which has exact chiral symmetry on the lattice. Con- 
sequently in contrast with the staggered Dirac opera- 
tor, the overlap operator has true zero modes even at 
finite lattice spacing. The non-chiral modes come in com- 
plex conjugate pairs. Again, we use the operator D°^D° 
which is Hermitian and positive definite and compute 
the lowest 5 eigenvalues in each chiral sector. The ma- 
trix sign function ejH^) is approximated by a 14th order 
Zolotarev expansion |25| 



(8) 



^ ' ^ I/tut 



with maximum errors ^ 0(10 ^") in the interval 
[0.04,1.5]. 



IV. SIMULATIONS 

Simulations are done with both the Wilson gluon ac- 
tion and the tadpole improved gauge field action for a 
range of lattice spacings (or couplings) and volumes. The 
physical scale is defined through the quenched string ten- 
sion -^/ct ~ OAAGeV. A new smearing method using uni- 
tarized Fat? links [26,] is used in the computation of the 
Wilson loops which substantially reduces the statistical 
fluctuations in the confining potential. Simulation pa- 
rameters and the measured string tensions a^cr are listed 
m TableHl The couplings were chosen carefully such that 
lattice spacings agree between simulations using the Wil- 
son gluon action and the improved gauge field action. A 
total of 1000 configurations were generated in each case. 
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FIG. 1: Spectral graphs for the different operators where chi- 
rality \x\ is plotted against the eigenvalue A. The ensem- 
ble is generated by the standard Wilson action at /? = 5.85 
(a ~ 0.123/m) on 10* lattice. Results are shown for 50 con- 
figurations. Note that the scale for A increases when the level 
of improvement increases. 
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FIG. 2: Distribution of the ratios of eigenvalues between the 
smaUest non-chiral mode and the largest "zero mode" for con- 
figurations with non-trivial topology. Results are shown for 
unimproved gauge fields at /? = 5.85 (a ~ 0.123/m) and lat- 
tice volume V = 10*. 

V. RESULTS 
A. Staggered Operator Improvement 

We first examine tlie effect of staggered operator im- 
provement on tlie infrared eigenvalue spectrum. Compar- 
isons are done with unimproved gauge fields for /? — 5.85 
(a ~ 0.123/m) and V = 10''. The same comparisons 
with improved gauge fields (same lattice spacing and vol- 
ume) will be presented in the next section. Results are 
shown in Fig. ^ where chirality |x| is plotted against 
the eigenvalue A for the different operators. Note that 
the chirality of an eigenmode is given by the expectation 
value of 75, which is a 4-link operator in the staggered 
basis . The same scale is used for the chirality. It can 
be observed that the eigenvalue spectrum depends quite 
sensitively on the way in which the staggered Dirac oper- 
ator is improved. Eigenmodes with small eigenvalues and 
relatively large chirality appear as the level of improve- 
ment increases. These small eigenmodes are taken to be 
the chiral "zero modes" which contribute to the index 
associated with the topology of the background gauge 
fields 1^ 0, . They can be identified by their separa- 
tion in energy and chirality from the rest of the spectrum 
(see also, Ref. [Gj). Throughout this project the follow- 
ing criteria are used for an eigenstate to be identified as 
a "zero mode" : i) it is at least a factor of two smaller 
in eigenvalue than the smallest non-chiral mode and ii) 
the chirality is at least five times larger than that of the 
smallest non-chiral mode. 

One can see also that as the level of improvement 
increases, the continuum fourfold degeneracy emerges 
where the scattered eigenmodes begin to form quartets. 
Note that at this lattice spacing with unimproved gauge 
fields the Asqtad operator is not sensitive to the topology 
and lattice artifacts are still dominant. One only starts to 








12 3 4 

IQI 

FIG. 3: Topological charge distributions (as a percentage of 
configurations) obtained by using different operators (unim- 
proved gauge fields, /3 = 5.85 (a ~ 0.123/m) and V = 10*). 

see a separation between the "zero modes" and the non- 
chiral modes with further improved operators. In addi- 
tion, a large renormalization is observed for the chirality 
of the would-be zero modes, which is |x| 0.5 instead 
of unity. This large renormalization was also observed in 
previous studies 0,113. Nevertheless, the "zero modes" 
can be identified without any difficulty in most of the 
cases for the UFat7 x Asq and HYP operators. 

To quantify the separation between the "zero modes" 
and the non-chiral modes, the distribution of the ratios 
of eigenvalues between the smallest non-chiral mode and 
the largest "zero mode" is plotted in Fig. [21for configura- 
tions with non-trivial topology (i.e., "zero modes" exist). 
We should mention that the configurations in Fig. |21are 
different for different operators as the topological indices 
obtained by using different operators do not always agree 
on a configuration by configuration basis (see next para- 
graph). It is possible that for a given configuration there 
are no chiral modes for one operator while they exist for 
the others. For comparison results of the overlap opera- 
tor are also shown. Theoretically, the ratio is infinite for 
the overlap operator because exact zero modes exist on 
the lattice for overlap fermions. It is finite here solely be- 
cause of computational precision. Results here show that 
the staggered operators are less sensitive to the topol- 
ogy at this lattice spacing with unimproved background 
gauge fields. The ratios are always three orders of mag- 
nitude for the overlap operator but the distribution is 
peaked at a ratio ~ 5 for the staggered Dirac operators. 
Nevertheless about 85% of the configurations have ratios 
> 5 which allows the "zero modes" to be identified. In 
addition our results also show that improvements using 
UFat7 links or hyper-cubic blocking are equally efficient. 

The identification of the topological charge index for 
a given gauge filed configuration is not unique. Because 
of lattice artifacts and the different ways the staggered 
operator is improved, it can be expected that the topolog- 
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FIG. 4: Same as Fig. 0but with improved gauge field action 
{P = 8.26 (a ~ 0.121/m), V = lO"). Note that the scale for 
A is larger in this case. 



ical charge indices obtained by using different operators 
do not always agree on a configuration by configuration 
basis. It is only when one is close to the continuum limit 
that the results agree. Even with the overlap operator 
previous studies showed that the topological indices 
do not agree on a configuration by configuration basis 
when, for example, a different mass parameter is used in 



Fat7XAsq — »- 

HYP "£> 
Overlap + 




Ratio 

FIG. 5: Same as Fig. |5|but with improved gauge field action 
{(3 = 8.26 (a ~ 0.121/m), V = 10*). 



the kernel. In the present case, we find that the topolog- 
ical charge indices determined by the different operators 
agree about 60-70% of the time, compared to 28% if the 
values were completely random. More explicitly, indices 
from different operators are compared pairwise and the 
agreement is 63% for overlap and UFat7 x Asq, 68% for 
overlap and HYP, 78% for UFat7 x Asq and HYP. Note 
that the agreement between UFat7 x Asq and HYP is 
larger than that obtained for an improved staggered op- 
erator and the overlap operator. It is then important 
to check whether the charge distributions are also dif- 
ferent because physical observables, e.g., the topological 
susceptibility, are related to the ensemble average of the 
topological charge. The distributions obtained by these 
operators are given in Fig. 0| It can be observed that 
there is no significant difference among the results. This 
is important because it indicates that physics is inde- 
pendent of the discretization scheme and one would ex- 
pect topological quantities obtained by these operators 
to agree. 



B. Improvement in the Gauge Field Action 

Lattice artifacts can be further suppressed if improve- 
ment is also applied to the gauge field action. This can 
be seen in Fig. ^where the infrared eigenvalue spectra of 
the different staggered operators are shown for configu- 
rations generated using the tadpole improved gauge field 
action. Here (3 = 8.26 (a ~ 0.121/m) and lattice volume 
is = 10^ so that both the lattice spacing and physical 
volume are very similar to those used in the unimproved 
case. Results here show that better topological properties 
are realized when the gauge field action is also improved. 
In particular even the Asqtad operator is sensitive to the 
topology at this coarse lattice spacing and "zero modes" 
can be identified unambiguously for the UFat7 x Asq and 
HYP operators. 

To have a quantitative picture of how improvement in 
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FIG. 6: Same as Fig. |21but with improved gauge field action 
{f3 = 8.26 (a ~ 0.121/m), V = lO'*). 



the gauge field action affects tlie spectral flow, we again 
plot the distribution of the ratios of eigenvalues between 
the smallest non-chiral mode and the largest "zero mode" 
in Fig. [SI for configurations with non-trivial topology. 
The same scale is used in Fig. |21and Fig. [S| In compari- 
son to Fig. 121 it can be observed that, for UFat7x Asq and 
HYP, the distribution has shifted significantly toward ra- 
tios of two or three orders of magnitude. Overall, more 
than 95% of the configurations have ratios > 5, com- 
pared to only 85% in the case with unimproved gauge 
fields. Hence gauge field improvement clearly increases 
the separation between the chiral "zero modes" and the 
non-chiral modes. We should again emphasize here that 
the ratios should be infinite theoretically for the overlap 
operator since exact zero modes exist on the lattice. 

We have also compared the topological indices ob- 
tained by the different operators in this case. It is found 
that the agreement increases significantly when the gauge 
field action is also improved: 91% between overlap and 
UFatTx Asq, 90% between overlap and HYP and 96% be- 
tween UFat7 X Asq and HYP which is again the highest. 
The charge distributions are given in Fig. |S1 It can be 
observed upon comparing with Fig. 13 that better agree- 
ment is obtained with improved gauge fields. Results 
here are significant because they indicate that different 
operators do respond the same way to the topology of 
the background gauge fields when discretization errors 
and lattice artifacts are reduced. In particular, results 
here show that the charge indices obtained by using the 
staggered operators and the overlap operator, two com- 
pletely different representations of the Dirac operator on 
the lattice, agree even on a configuration by configuration 
basis at a high percentage as the level of improvement in- 
creases. 



C. Dependence on Lattice Spacing 

Here we examine the dependence of the infrared eigen- 
value spectrum on lattice spacing and study the spectral 
flow as one approaches the continuum limit. Calculations 
are done at three lattice spacings with fixed physical vol- 
ume V ^ (1.2/m)^ (so lattice volume increases as lattice 
spacing decreases, see Table IJ). Results are shown for the 
Asqtad and UFat7 x Asq operators in Fig. 7. One sees 
that separation between the "zero modes" and non-chiral 
modes becomes more clear. In addition the continuum 4- 
fold degeneracy is better realized as one approaches the 
continuum limit. Note that the chirality of the "zero 
modes" increases as the lattice spacing decreases and it 
is larger for the UFat7 x Asq operator. This gives ev- 
idence to the fact that discretization errors and lattice 
artifacts are indeed responsible for the failure of stag- 
gered fermions to be sensitive to gauge field topology on 
coarse lattices. The chiral "zero modes" associated with 
the topology of the background gauge fields emerge as 
one approaches the continuum limit, a conclusion also 
obtained in 

Our results also show that it is necessary to use a lat- 
tice spacing a < 0.1/m for the Asqtad operator to be 
sensitive to gauge field topology. This is similar to the 
lattice spacings used nowadays in dynamical simulations 
of QCD using the Asqtad fermion action and the tadpole 
improved gauge field action On the other hand, 

further improvement of the staggered quark action, e.g, 
actions with unitarized fat-links UFat7 x Asq, may be 
necessary when working on coarse lattices to ensure that 
continuum physics is reproduced correctly. 

D. Dependence on Lattice Volume 

As discussed in the introduction, we expect a strong 
dependence of the spectral density on the volume V of 
the lattice and the topological charge average (Q^) should 
scale with V. In Fig. 8 the infrared eigenvalue spec- 
trum of the UFat7 x Asq operator is given for lattice vol- 
umes V — 8'^, 10^, 12^ and 16* at a fixed lattice spacing 
a — 0.121/m (/3 = 8.26) with improved gauge fields. Re- 
sults from previous sections showed that "zero modes" 
should be visibly separated from the non-chiral modes 
at this spacing for configurations generated by the im- 
proved gauge action. This is the case when the volume is 
~ (1.2/771)"*. A similar volume effect has been seen with 
the overlap fermion operator also 0, ■ A reasonable 
assumption is that this minimum volume effect is related 
to a property of the gauge field, perhaps some minimum 
volume is necessary before the topological structure is 
fully formed. 

As the lattice volume increases the number of eigen- 
values increases and one may expect the gap between the 
zero modes and the non-chiral modes to scale as [2^ . 
The net result is that there are more low-lying non-chiral 
modes as volume gets bigger. As Fig. 8 shows, a certain 
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FIG. 7: Dependence on lattice spacing. Results are shown for the Asqtad (top) and the UFat7 x Asq (bottom) operators with 
improved gauge fields. Physical volumes are ~ (f.2/m)* in all cases. 
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FIG. 8: Spectral graphs for the UFat7 x Asq operator at a fixed coupling f3 = 8.26 (a ~ 0.121/m) for volumes V = 8**, 10'', 12** 
and 16* with improved gauge fields. The spectrum becomes noisy again for large volumes. Note that the scale for A decreases 
with increasing volume. 



volume is necessary before the would-be "zero modes" 
show up and sensitivity to topology is established. As 
the volume is increased even further the number of chiral 
"zero modes" increases roughly as Vv as can be inferred 
from Fig. |51 In the configuration average, the tail of the 
chiral-mode distribution begins to merge with the low 
eigenvalue tail of the non-chiral mode distribution. How- 



ever, the indentification of the would-be "zero modes" is 
not quite as difficult as it may appear from Fig. 8. As 
seen from Fig. IHlthe criteria which we adopt for a "zero 
mode", namely, at least a factor of 2 smaller in eigen- 
value than the smallest non-chiral mode {Rx = 2) and a 
factor of 5 larger in chirality (i?^ = 5), are very robust 
even for the IG** lattice. Imposing other values for these 
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FIG. 9: The square-root of the total number of zero modes 
of the UFat7xAsq operator as a function of lattice size L 
(lattice volume V — L'*). Results of using different criteria 
for identifying the zero modes are shown. 




FIG. 10: Topological charge distribution (as a percentage 
of configurations) for different lattice volumes. Results are 
shown for the UFat7 x Asq operator at l3 — 8.26 (a ~ 
0.121/m) with improved gauge fields. 



factors gives very similar results. 

The expected increase of (|Q|) with the volume of the 
lattice can be seen in Fig. ^|where the topological charge 
distribution is shown for different lattice volumes. Most 
configurations have trivial topological structure when the 
volume is small and the charge average {\Q\) increases 
gradually with the volume of the lattice. 



VI. COMPARISON WITH RMT 

Results from our systematic analysis give strong in- 
dications that staggered fermions are sensitive to gauge 
field topology when lattice artifacts are sufficiently sup- 
pressed. It is then important to compare the distribution 
of the non-chiral modes, after successfully identifying the 



would-be zero modes, against the predictions of RMT. In 
particular, we study the cumulative distribution of the 
smallest non-chiral modes in different topological sectors 
PcSlmiC) where RMT gives the predictions |5l.l30| 



xe 



]: (\h{xf ~ h{x)h{x)\ e'^'/*dx.{9) 
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Here, P'^^^x) denotes the distribution of the smallest 
non-chiral modes in charge sector Q and In (a;) is the 
modified Bessel function of the first kind (order n). The 
variable ^ = AEF is related to the infinite-volume chiral 
condensate E which is a fitting parameter. 

Calculations are done for three different operators 
((UFat7) X Asq, HYP and overlap) on two ensembles 
of configurations (unimproved and improved gauge fields 
with a ^ 0.12/m, V = IC*). These ensembles are specif- 
ically chosen because the size of the lattice is limited 
by the cost in computing the eigenvalue spectra of the 
overlap operator while a coarse grid is used so that the 
physical volume is larger than the critical value (1.2/m)*. 
This criterion is important as it ensures that we are in 
the e-regime where RMT is applicable ^llj . 

In our analysis, the would-be zero modes are first iden- 
tified and the configurations are classified according to 
their charge indices obtained using the index theorem 
(Q. We should re-emphasize here that the indices ob- 
tained by using different operators do not always agree 
on a configuration by configuration basis but the charge 
distributions are indistinguishable (see Fig. ISandEJl. It 
should also be noted that all configurations are used in 
our analysis. This is different from some previous stud- 
ies 1^ @ where a small portion of the ensemble was dis- 
carded. In these studies the topological charge indices 
were also calculated using a direct discretization of the 
continuum formula, Qg ~ / ^^"^'^ F^^Fpcr- Only those 
configurations with approximate integer value of Qg were 
included in the analysis in Ref. ^ , and configurations for 
which the charge index obtained by the index theorem is 
different from Qg were excluded in Ref. |^. 

Results are shown in Fig. ^2 The solid curves are 
predictions from RMT. Because computing the spectra 
for the overlap operator is much more expensive, the 
sizes of the ensembles arc different: 1000 configurations 
for (UFatT) x Asq and HYP but only 400 for the overlap 
operator. This is the reason that the overlap results 
appear to have poorer agreement with the predictions of 
RMT. Note that the agreement is better with improved 
gauge fields. The results given here contribute to the 
evidence that staggered fermions do feel gauge field 
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FIG. 11: Cumulative distribution of the smallest non-chiral modes in different topological sectors. Results are shown for 
unimproved (top) and improved (bottom) gauge fields at a ~ 0.12/m, V = lO"*. Solid curves are predictions of RMT. 



topology, provided that lattice artifacts are suppressed 
considerable, as the agreement with the predictions of 
RMT is impressive. They should be compared with 
previous studies 0, with the unimproved staggered 
operator which appeared to indicate the presence of only 
a trivial topological sector. 



VII. CONCLUSION 

In this project we studied numerically the spectral 
properties of a variety of improved staggered Dirac op- 
erators. Four systematics have been examined: i) im- 
provement in the staggered operator, ii) improvement in 
the gauge field action, iii) lattice spacing and iv) lattice 
volume. 

It has been observed that the infrared eigenvalue spec- 
trum depends sensitively on the way in which the stag- 
gered fermion operator is improved. On coarse lattices 
the unimproved operator is insensitive to gauge field 
topology. As the level of improvement increases, either on 
the operator itself or the background gauge fields, eigen- 
modes with small eigenvalues and large chirality appear. 
These small eigenmodes can be identified as the chiral 
"zero modes" associated with the topology of the gauge 
fields. Sensitivity to the topology also increases as one 
approaches the continuum limit. This gives evidence that 
lattice artifacts are responsible for the failure of the unim- 
proved staggered operator to reflect properly the gauge 
field topology on coarse lattices. Our results also show 
that a lattice spacing a < 0.1/m is enough for the Asqtad 



operator to have a correct response to the topology with 
improved gauge fields. This spacing is of the order of the 
lattice spacings used in present day state-of-the-act dy- 
namical simulations of QCD. On the other hand, the next 
level of improved staggered operators, e.g., UFat7 x Asq, 
may be required to produce configurations which describe 
the correct continuum physics on coarser lattices. We 
also observe that the topological charge distribution is 
independent of which operator is used even though the 
charge indices do not always agree on a configuration by 
configuration basis. However, the agreement increases 
with the level of improvement. A minimum physical vol- 
ume of about (1.2/m)'' seems to be necessary in order for 
"zero modes" to show up and for sensitivity to topology 
to be established. This effect observed earlier 

for overlap fermions. As volume is increased the number 
of would-be "zero modes" increases as ^/V but the num- 
ber of low-lying non-chiral modes increases faster so some 
merging of tails of the distributions takes place. Up to 
the 16^ volume considered here it is still possible to make 
a clear separation of would-be chiral modes from non- 
chiral modes. The distribution of the non-chiral modes 
is matched with the predictions of RMT. The agreement 
is comparable to that obtained using overlap fermions. 



Based on this, and other work done in the past 

year, one has strong evidence that, provided one uses 
improved staggered operators and improved gauge fields, 
staggered fermions properly feel gauge field topology. 
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